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Abstract 

We construct a Dark Matter (DM) annihilation module that can encompass the predic- 
tions from a wide array of models built to explain the recently reported PAMELA and 
ATIC/PPB-BETS excesses. We present a detailed analysis of the injection spectrums 
for DM annihilation and quantitatively demonstrate effects that have previously not 
been included from the particle physics perspective. With this module we demonstrate 
the parameter space that can account for the aforementioned excesses and be compat- 
ible with existing high energy gamma ray and neutrino experiments. However, we find 
that it is relatively generic to have some tension between the results of the HESS exper- 
iment and the ATIC/PPB-BETS experiments within the context of annihilating DM. 
Wc discuss ways to alleviate this tension and how upcoming experiments will be able 
to differentiate amongst the various possible explanations of the purported excesses. 



1 Introduction 



Recently there has been a series of experimental results suggesting that we may have indi- 
rectly detected dark matter (DM) within our Galaxy. The combination of the positron frac- 
tion measured by the PAMELA experiment [1] and the ATIC/PPB-BETS experiments [2,3], 
have led to a compelling picture of DM being responsible for a new population of positrons 
at high energies. These excesses, if confirmed, could in principle have alternative explana- 
tions through either refining our understanding of charged particle propagation within our 
Galaxy, or by identifying new astrophysical sources of positrons coming, for instance, from 
pulsars [4]. It is intriguing therefore, that new experimental data expected in the near future 
could not only confirm or contradict those results, but also allow us to possibly determine 
the physics behind these excesses. Furthermore, in the case of DM, such experiments could 
strongly constrain the various DM models. 

Broadly speaking, the plethora of DM models bifurcate into either annihilating [5, 6] or 
decaying [7] DM. In this paper we choose to focus on the former possibility as being the 
source of the electronic excesses. The above experiments then place strong restrictions on 
the models, so we adopt the following phenomenological inputs as constraints: 

• There is an excess in the flux ratio $e+/^e++e- observed by the PAMELA experiment 
extending to at least 100 GeV [1]. 

• There is an excess in the ATIC/PPB-BETS experiments for the flux of charged elec- 
trons and positrons, $e++e-, extending to energies of ~ 700 GeV [2,3]. 

• There is no excess observed by the PAMELA experiment in the antiproton flux [8] . 

• In the absence of large local overdensity in the DM distribution (boost factor), the 
annihilation cross section in our galaxy needs to be O(IOO) times larger than a standard 
thermal WIMP. 

These facts arc not easily reconciled. The last of these assumptions follows from the 
large measured rates combined with the higher mass scale indicated by the ATIC/PPB- 
BETS anomaly. For the case of a WIMP DM, a large enhancement of the cross section is 
needed [9-11]. Alternatively, a large boost factor (BF) is required. However, such a pos- 
sibility seems unlikely in light of the results from N-body simulations [12]. Furthermore, a 
model must prefer annihilation into leptonic final states so that the antiproton fraction is 
not overpopulated. 

There has been a recent explosion in model building that attempts to incorporate the 
necessary ingredients to explain these excesses. Typically, these models explain the electronic 
activity by either assuming a symmetry that forbids hadronic production, or otherwise pos- 
tulating an intermediate light state that can only decay into light leptons due to kinematics. 
Most of these studies have either stopped at the heuristic level of explanation, or attempted 
quantitatively only to postdict certain experiments. It is therefore desirable to consider 
a larger set of experimental data in order to better establish the correct model-building 
direction. 
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We attempt to address the following questions: 

• Given a model that can explain the PAMELA and ATIC/PPB-BETS data, what are 
the experimental bounds arising from other experiments? 

• What are the viable classes of models? 

• For these models, what are the implications for upcoming experiments? 

The most logical additional signature which has not been entirely explored is the one 
coming from photons. Whenever there are charged particles in the final state there will 
be additional photons radiated, leading to a model independent signature [13]. Additional 
sources of photons may contribute depending on the specific details of the model. Recently 
there have been a few papers [14-16] that have studied the bounds from high energy photons 
in models that explain the excesses. The authors of [15, 16] reached the conclusion that for 
most dark matter density profiles, experimental results rule out the possibility of annihilating 
DM as an explanation of the excesses. In [14,15] the case where DM directly annihilates into 
a pair of SM leptons was studied. We reach a similar conclusion to [15], that such models 
disagree with the experimental data collected for high energy photons. In [16] models where 
DM annihilates through a light state and then into leptons was studied. In our paper we 
focus on these models and reach a different conclusion than the authors of [16]. While we find 
there exists some tension between models that explain ATIC/PPB-BETS and high energy 
photons, they are not ruled out by an order of magnitude. We also include several effects 
that have not yet been studied in the hterature, that can amehorate this tension, such as 
dark sector radiation. 

To study the implications of the present experimental data, we construct a module that 
incorporates many of the required features necessary for a model to explain the excesses. 
This module has several parameters that allow us to interpolate between different classes of 
models. 

Our main results are summarized as follows. For annihilating DM scenarios that explain 
both PAMELA and ATIC/PPB-BETS: 

• Rather generically, such models are in tension with constraints from high-energy pho- 
tons. 

• Photons are more constraining than the antiprotons measurements. In particular, 
models that produce antiprotons and still fit the PAMELA data in many cases produce 
too many photons to be consistent with the measurements. 

• The tension is not sufficient to exclude all models. However, it requires a factor of 
order 0{2 — 5) that may arise from various sources, e.g. a local boost factor or a less 
cuspy DM profile. 

We find a number of interesting implications based on these results. Upcoming experiments 
have the power to exclude the full region of parameter space, for models that explain the 
ATIC/PPB-BETS excesses with annihilating DM. From this point of view decaying DM 



2 



models are attractive because they generate fewer photons at the center of the galaxy. Models 
of annihilating DM that do not seek to explain the ATIC/PPB-BETS are also viable, and 
can be tested with several ongoing experiments. In particular, we stress the importance of 
the currently running experiment, FERMI [17], for helping determine the underlying nature 
of these excesses. 

This paper is organized as follows. In Section [2] we construct a DM module that enables 
us to investigate the parameter space of models. We then calculate the particle physics 
input for the relevant experiments, namely the injection spectrums for e"*", p, 7 and p. In 
Section [3] we discuss the astrophysical inputs for our study. We review the methods used in 
this study for propagating the various particles from their source to Earth. Additionally, we 
further discuss the experimental inputs that we use, and comment on many of the uncer- 
tainties associated in calculating the fluxes for them. In Section |4] we present the resulting 
astrophysical fluxes calculated from our particle physics module. We demonstrate how the 
various particle physics and astrophysics parameters affect the predicted fluxes for e"*", p, 7 
and u. Finally, in Section |5] we discuss the implications and interpretations of the regions of 
parameter space that we find to be consistent with the experiments studied. In Appendix [A] 
we discuss the so called "leaky box" approximation which we use to estimate the positron 
background flux. 

2 Unified Dark Matter Module 

We are interested in understanding the bounds and predictions for characteristic DM models 
that could be used to explain the excesses observed. In this work we do not focus on 
the bounds for a particular model nor are we completely model independent. Instead, we 
construct a module that contains the most important components that we identify from the 
particle physics perspective. This module can then be appropriately recast to reflect the 
predictions from a wide array of models that have been, and inevitably will be, built. 

We construct the dark matter module loosely in accord with the Unified Dark Matter 
model of Arkani-Hamed et al. [5]. This model offers the intriguing possibility to describe 
not only the PAMELA and ATIC/PPB-BETS excesses, but also the INTEGRAL [18] and 
DAMA [19] excesses simultaneously. We treat [5] as representative of a class of ideas, and 
choose to incorporate those features of [5] that are relevant to studying the experimental 
consequences for indirect searches in the high energy e"*", p and 7 channels. 

The components we choose to include in our module for explaining the excesses are 
the following. We assume there are heavy DM particle(s), x(s), that are charged under 
some "dark" gauge group, and possibly the SM electroweak (EW) gauge groups as well. 
Additionally, the dark gauge group is broken and therefore the sector consists of light gauge 
bosons which we collectively refer to as 0. The light gauge bosons are required for two 
reasons: On the one hand, they allow for a kinematical explanation for the electron but no 
antiproton excess measured by the above experiments. On the other hand, they play part 
in the Sommerfeld mechanism that can enhance the usual thermal WIMP cross-section to 
the required rate today. With this in mind, we allow the x's to annihilate either through 
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the SM EW gauge bosons, V = {W, Z), or the 0's: 



XX yy or XX 00- 



(2.1) 



If X annihilates through SM gauge bosons, then the final states are clear. The annihilation 
into needs further explanation. If one assumes that there are measurable consequences in 
experiments then ultimately needs to decay into SM final states. There are two possibilities 
that allow to decay. Either the SM matter fields are charged under the dark gauge group, 
in which case can decay into SM fields based on their charge assignments. Otherwise, the 
gauge bosons can mix with the SM gauge bosons and thereby decay through this mixing. 
While either scenario is possible in principle, we will choose the latter and couple to the 
SM matter through gauge boson mixing. On general grounds the vector will decay back 
to SM states by mixing through the photon or the Z boson. However, given the lightness of 
0, the decays going through Z mixing will be further suppressed by at least m^/m|. Unless 
the 7-0 mixing is much smaller than the Z-0 mixing, one can then assume that couples to 
SM particles proportionally to their electric charge. Hence the decay branching fractions 
are completely determined by its mass. An example annihilation is shown in Figure [TJ 



Figure 1: An example of x annihilation, xx — ^ 00 — > e'^e e^e via the mixing of and 7. 

With these considerations we have come up with a minimal module that has five particle 
physics parameters: 

• - sets the mass scale for the DM annihilation. 

• - determines what states can decay into and the kinematics of the decay products. 

• «DM - the strength of the gauge coupling in the dark sector. 

• {av) = {<7v)vv + {'^'^)<i>^ - ^ free parameter for the overall cross section. 

• RsM = ^T^'C'^ - a free parameter for the relative contributions of the annihilations into 





the SM vs dark gauge bosons. 
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While some of these parameters may seem redundant, they allow us to cover the parameter 
space of a large number of models without having to calculate within each model separately. 
In particular, the inclusion of odm as a separate parameter is noticeable as we keep {av) and 
RsM free. 

Once we stipulate that is a dark gauge boson, we need to allow for the possibility 
that the dark gauge group is nonabelian. In fact this is exactly what is desired in [5] to 
explain DAMA and INTEGRAL anomalies. If represents collectively the gauge bosons 
of a nonabelian group, then there are additional processes for annihilation into the SM 
compared to those shown in Figure [ij Just like in QCD once the gauge bosons are produced. 




Figure 2: An example of XX annihilation, with the inclusion of the dark sector shower. 

they can shower and split into new gauge bosons. We give an example of this process in 
Figure [2} where after showering the dark gauge boson subsequently decays into SM final 
states through photon mixing. This process has not been quantitatively explored before in 
this context and we demonstrate the effect in later sections. As we shall see, while naively 
there is no large enhancement as the 0's are massive, a significant change in the resulting 
energy spectrum arises since we assume ^ m^. 

There additionally could be another annihilation channel for x, i.e. xX ~^ ^'P- This is 
strongly dependent on the model, and could in principle be used as a separate parameter in 
our module. For instance this mode will not occur with any appreciable rate in models where 
X is a single Majorana particle. We choose not to include this as a separate parameter, and 
instead one can infer bounds on this mode from our Rsm appropriately rescaled. 

2.1 Calculation of Injection Spectrum 

In this section we calculate the particle physics input for all dark matter indirect experiments. 
For the experiments that we are interested in, we simply need the injection spectrum of e"*",/?, 
7, and v coming from the annihilation of the x's. Therefore we need to calculate the inclusive 
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annihilation of XX threshold and then extract the energy distribution, dN/dE, for each 
particle we are interested in. 

To calculate dN/dE we implement our module in several MC programs and scan over 
the particle physics parameters m^^m^ and odm that define it. We do not calculate (erf) 
and RsM from first principles, instead we will fit the experimental data for PAMELA and 
ATIC/PPB-BETS to fix these parameters. The result of this fit will then dictate what 
boost factor is needed compared to a standard thermal WIMP annihilation cross section. If 
one calculated the Sommerfeld effect within a model, this would constrain the values of the 
masses and gauge couplings given the necessary cross section. However, since we do not wish 
to focus on one particular model alone, we do not require this consistency check. Additionally 
once the mass and odm are given, a model could predict the ratio Rsm- However, Rsm can 
also be an independent of odm, for instance if the DM fields x are singlets under the SM 
and mixed through a Yukawa coupling. 

To implement the module in event generators we make the specific choice of SU (2) for 
the dark gauge group, and we define x to be a bi-doublet under the SM SU{2) and the 
hidden SU{2). This choice only affects the dark gauge boson parton shower directly, since 
(av) and Rsm are fit to data. However, since we will scan over odm, a different choice of 
coupling can still be be used to approximately interpolate amongst different Casimirs and 
extrapolate the results to different gauge groups. 

The parameter space we choose to cover is shown in Table [T} Our choice of whether to 
implement a parton shower in the dark sector is a binary one. This allows us to cover the 
case where the dark gauge group is f/(l) and the gauge bosons don't shower. 



No dark sector parton shower 






200 GeV - 3.5 TeV 
200 MeV, 500 MeV, 1.2 GeV, 8 GeV, 15 GeV 
arbitrary 


Dark sector parton shower 




ttDM 


200 GeV - 3.5 TeV (discrete grid) 
200 MeV, 500 MeV, 1.2 GeV, 8 GeV, 15 GeV 
10-^ 10-2, 4- 10-2, 10^1 



Table 1: Parameters scanned in the DM module. 

To implement the decays of the particle we use different effective field theories depend- 
ing on the mass of (p that we are interested in. When ^ ^qcd, we coupled (j) directly 
to the quarks at parton level, which can then be showered and hadronized. On the other 
hand if < a few GeV, this procedure is a bad approximation. In this case we couple 
(j) directly to mesons. We take the corresponding cross sections from experimental data of 
e~^e~ — i> hadrons exclusive processes [20]. In particular for the 500 MeV case, we assume 
decays to e^e~, ^tt'^tc' and the ratio is approximately 2:2:1 as shown in Table [2j For 
the 1.2 GeV case we include other mesonic 2-body modes as well as 3- and 4-body decays 
that are non- negligible. We did not implement any other intermediate mass between 1.2 GeV 
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and 8 GeV because multi-body final states become increasingly important and there is not 
enough experimental information to completely reconstruct the exclusive final states. The 
8 GeV and 15 GeV were chosen to have separated enough from the quarkonia resonance 
region (where the hadronization model we use from Pythia [21] suffers from large uncertain- 
ties) and be above and below the bb threshold. We catalogue the decay modes implemented 
and branching fractions of in Table [2] 











Mode 


BF 




Mode 


BF 




8 GeV 


0+ — » e+e- 


1.6 - 10-1 


200 MeV 


e+e^ 


1 






^ fJ'^fJ'^ 


1.6- 10-1 


500 MeV 




4- 10-1 






— > r^r 


1.6 - 10-1 




1 _|_ 
yU^/i 


4- 10-1 






— > MM 


2.1 - 10-1 




— > TT+TT- 


2 ■ 10-1 






(f) —>■ dd 


5.2 - 10-2 


1.2 GeV 


e+e^ 


3.4- 10-1 






(f) ^ cc 


2.1 - 10-1 




fi^jJ,^ 


3.3-10-1 






— ^ ss 


5.2 - 10-2 




— i> UJTV^ 


7.9 • 10-2 




15 GeV 


0+ ^ e+e- 


1.5- 10-1 




— ^ 71^71-71*^71*^ 


7.5-10-2 






— i> jl^ jl^ 


1.5- 10-1 




— > 7r+7r- 


6.4- 10-2 






— s> r+r- 


1.5-10-1 




^ K+K- 


4.5 - 10-2 






(j) uu 


2.0 - 10-1 




— > 7r+7r+7r-7r- 


4.1 - 10-2 






(f) ^ dd 


5.0-10-2 




— > 7r^7r-7r'^ 


2.4- 10-2 






(p ^ cc 


2.0 - 10-1 




^ 


5-10-3 






(f) ^ ss 


5.0- 10-2 










(j)^bb 


4.8 - 10-2 



Table 2: Branching Fractions for 0. The values for 0.5 and 1.2 GeV are extracted from 
experimental data for exclusive e'^ + e~ —>■ hadrons processes [20] , while for 8 and 15 GeV 
are computed using BRIDGE. 

To calculate the injection spectrums using existing Monte Carlo (MC) tools is quite dif- 
ficult. The kinematic regime we study is based on very heavy particles annihilating through 
very light particles, that subsequently decay. In this regime most MC generators that we 
have used have difficulties. To generate our injection spectrums we were forced to use a 
variety of generators linked together depending on the task: MadGraph/MadEvent [22], 
BRIDGE [23], SHERPA [24], and Pythia [21]. 

For the SM annihilations XX ~^ ~^ fermions, we generate parton level events keeping 
spin correlations. We then shower and hadronize, including the effects of photons that are 
showered from the W gauge bosons [25]. To generate xX — > 00 — * SM without the dark 
sector parton shower we again generate parton level events including spin correlations. We 
then shower and hadronize either the fundamental particles, or just shower the charged 
mesons and their decay products. Unlike [16], when applicable, we include the C(l) effects 
of calculating photons showered from muon decays. 

For the case when we include the dark sector parton shower, we first need to calculate the 
massive vector boson splitting function for the SU{2) case that we have implemented. Given 
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the kinematics that we scan over, ^ m^, the effects of the massive splitting function are 
minor and it is sufficient to use the massless splitting functions, 



aoM^il - z{l - z)y 



2n 



z(l-z) 



(2.2) 



Nevertheless in our computation we include the complete massive splitting function which 
only differs by the inclusion of another term that is subdominant over most of our kinematic 
range. To generate events including the dark sector parton shower, we compute the 2 — 2 
matrix element. We then shower in the dark sector and decay the keeping some of the 
spin correlations. Finally, we shower and hadronize the SM particles that come from the 
decays, including the effects listed previously. 



Ainiproloii Speclni 




0.005 0.010 



0.050 0.100 



0500 1.000 




0.005 0.010 



0.050 0.100 



0.500 1.000 ' 



ojoo I mi 




0.005 0.010 



0.050 0.100 



OJOO 1.000 



Figure 3: dN/dx for p, e"'",7, and v for the various m<^ used. The SM contribution from 
XX ~^ W is also shown. In the range of x plotted we never reach Aqcd/^x ^^^s the shapes 
are universal and independent of m^, even for the antiprotons. 

In Figure [3] we plot dN/dx, where x = E/m^, for p, e+, 7, and v for different m<^ with no 
parton shower included. We additionally include the SM annihilation channel XX ~^ VV. For 
later reference, we also plot in Figure |4] the photon and positron spectrum for the two body 
decay xX ~^ e^e~ . Both the SM annihilation channels will be used below for comparison 
with different models and to study the constraints on annihilating though the SM gauge 
bosons. As one can see, the larger the number of open decay channels for 0, the larger the 
population of low energy particles, namely the softer the spectrum. One can see that in 



8 



Positron Spectrum, xx^^^^ ' ^;t'-*'^0->2f^2f? 

dNIdx 
10.00 r 

5.00 - 



1.00 
0.50 

0.10 
0.05 

-J O./O 0.15 020 030 050 0.70 /.OO ^ O./O 0./5 0.20 OJO ^ (L50 0.70 IW 

Figure 4: Photon and positron spectrum for the two body decay XX ~^ e+e^ (dashed red 
hne). For comparison, we also plot the four body decay xX ^ 20 ^ 2e^2e~ (solid blue 
line) . 

the 8 GeV, 15 GeV and SM gauge boson cases, showering and hadronization effects from 
QCD induce several orders of magnitude increase in dN/dx. Moreover, showering effects 
also softens the spectrum of the direct decay into e~^e~, as is clear from Figure |4j 

In Figure [5] we plot dN/dx, for p, e"'",7, and u for different values of odm to illustrate 
the effects of the dark sector parton shower. Without showering, the spectrum, dN/dx, is 
universal and independent of m^, this is no longer true when showering is taken into account. 
The reason for this is that the number of energy decades for showering depends on the ratio 
m^/m^. Hence in order to demonstrate the odm dependence, we fix the value of and 
m^. One should note that showering can have a large effect on the dN/dx. Hidden sector 
showering not only softens the spectrum via radiation, but also through additional decays. 

2.1.1 Particle Physics Uncertainties 

In this section we review some of the uncertainties in calculating the dN/dE for e^, p, v 
and 7. Most channels are rather clean from the particle physics perspective, but there are 
several issues that can have a significant effect on the dN/dx that we compute. These issues 
are hadronic uncertainties in the calculation of the dN/dx for p, showering uncertainties 
which change the photon spectrum and the effects associated with the possibility of having 
different matrix elements which can affect all particles. 

Hadronic uncertainties are due to the fact that once we have SM quarks we need to turn 
these partons into mesons and baryons through hadronization/fragmentation. While QCD 
can correctly describe the parton showering evolution, the process of hadronization relies 
upon phenomenological models that are tuned to data. While the effects of hadronization are 
important in some specific collider searches, when discussing high energy jets at colliders, one 
rarely talks about exclusive channels that label specific numbers of mesons and baryons. On 
the other hand, in the case of DM we look at the fully exclusive p channel. To hadronize we 
use the Pythia program. Pythia employs several possible fragmentation models for creating 
baryons such as diquark, popcorn and advanced popcorn. The Pythia manual [21], states 



Photon Spectrum. XX ^ • XX~^4'4''^^^ 
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Positron Specf rum. = 1 TeV. ^ 200 MeV 




Photon Spectrum, m ^ = J TeV. = 200 MeV 



0.005 0.010 



0.050 0.100 



Neutrino Spectrum, m^ = } TeV, m^ = 500 MeV 
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oms 0.010 
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0.005 0.010 
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0.500 1.000 



Figure 5: The effect of varying odm over tlie range given in Table [T] We plot dN/dx for a 
given particle type with a fixed niy. = 1 TeV and various m^. 



that the tuning of their baryon production models are done with a global fit to the data and 
that the resulting fragmentation functions for individual baryons can be lower or higher than 
the actual data. To quantify this effect we used the Pythia default popcorn algorithm, and 
compare it with the advanced popcorn algorithm and SHERPA's AHADIC++ algorithm. 
We find that there is an 0{1) effect in using different hadronizat ion/fragmentation algorithms 
for the p injection spectrum. 

For showering, there are various levels of sophistication that one may try to employ. For 
instance in comparison with [16] we include the effects of showering off the decay products 
of the muon while they only include the shower from the muon itself. We find that this is an 
0{1) effect over the entire range of energies that we examine. Additionally there are further 
uncertainties depending on how the parton shower is implemented. For instance, when 
SHERPA is used to compute the radiative decays of the vector 0, taking into account the 
full matrix element structure for /^/~7, it tends to underpopulate the high energy gammas 
compared to the parton shower approximation. Since the above effect is model dependent, 
we do not consider it here and instead, rely solely on the parton shower. Nevertheless, one 
should keep in mind that this discrepancy induces another 0{1) effect that may reduce or 
enhance our results. 

As in the case of decays, differences between matrix elements can introduce a nontrivial 
dN/dx shape dependence which translates into different fluxes observed in experiments. This 
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has been studied in the hterature e.g. [26] . Including non trivial effects of spin correlation, 
typically do not introduce large changes in the positron spectrum [26], since usually the 
spectrum changes only at very high energy. However, depending on the form of the positron 
injection spectrum, this can lead to non trivial changes in the region of interest for gamma ray 
searches. While we only take the matrix element for 2 — > 4 assuming there is an intermediate 
light gauge boson, one should keep in mind that there could be a factor of a few uncertainty 
if the underlying model was different. 



In this section we review the astrophysical propagation of e~^, p, v and 7. Due to a lack 
of theoretical understanding and experimental data, propagation models suffer from large 
uncertainties and there are various models that exist in the literature. The choice of a model 
is not only important for estimating the DM signal, but also for correctly evaluating the 
backgrounds. Using different models may result in significantly different predictions. 

As an example, we could, in principle, use a program such as GALPROP [27] and cal- 
culate both the background and signal propagation by using a best fit propagation model. 
However, we aim to be as data driven as possible and introduce the minimal amount of the- 
ory necessary to estimate both the signal and background. We therefore use semi-analytical 
propagation method which we review below. This approach has been studied extensively 
in the literature [28, 29] and it therefore allows for a simple estimation of the underlying 
uncertainties. In later sections we will briefly revisit some of these uncertainties. 

We can also employ a data-driven analytic approach to estimate the backgrounds. When- 
ever wc can, wc use known measurements or otherwise conservative models for the evaluation 
of the backgrounds. Below we discuss the propagation model, the experimental data and the 
background estimation for each of the relevant channels. Wc discuss in each case the uncer- 
tainties involved and explain how we take those into account when fitting the predictions to 
the data. 



In this section we review the DM density profiles that we use as inputs when calculating the 
fiuxes from DM annihilation. Most of the dark matter profiles that we consider here are in- 
ferred from N-body simulations. Starting in the mid-nineties, a paradigm emerged [30] where 
by examining the results of dark matter N-body simulations for many different galaxies, a 
type of universality for the density profiles appeared. This led to the famous NFW profile [30] 
for dark matter that is in common use today. Since then other groups have examined this 
universality and found similar results. Nevertheless, there is some disagreement between 
groups concerning how cuspy the dark matter profile is at the center of the galaxy. Most 
standard dark matter profiles can be parameterized using the {a, P, 7) parametrization [31] 



3 Astrophysics Inputs 



3.1 DM Halos 




(3.1) 
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(1) 



Addi- 



The two most commonly used profiles of this type are NFW [30] and Moore [32 
tionally the Isothermal profile [34] also fits into this parameterization, but it is based on older 
approximations for dark matter halos and there is no evidence for it in N-body simulations. 
While we do not consider the Isothermal profile viable, we show it below for the sake of 
comparison. In (3.2) we hst the values of {a,P,'y,rs) that correspond to the commonly used 
profiles. 



ProfileName 


a 




7 




NFW 


1 


3 


1 


20 


Moore 


1 


3 


1.16 


30 


Isothermal 


2 


2 





5 



(3.2) 



Recently N-body simulations have been able to increase their resolution by using O{10^) 
particle simulations. From this, one is able to test the validity of the commonly used profiles 
such as NFW. It turns out that, as the resolution has increased and one has been able to 
probe further and further in towards the center of the Galaxy, the cuspiness does not persist. 



This has led [35] to propose that, instead of the commonly used profiles in (3.2), the best fit 
profile is the so called Einasto profile [36 



(2) 



p(r) 



PQ exp 



a 



(3.3) 



What determines the deviation from the power law behavior is the Einasto a parameter 
which was found in [35] to have a mean value of 0.172 when fitting to simulations based on a 
few hundred million particles carried out by the Virgo consortium. Recently billion particle 
simulations have been carried out. The Aquarius [38], Via Lactea [39], and GHALO [40] 
simulations all fit the Einasto profile and a has been found to lie in the range [0.142, 0.177]. 
To be conservative, we choose a larger range of a consistent with [35]. Below, we will examine 
a = 0.12,0.17,0.20. 

The cuspiness of the profile is not the only important feature, the local density is obviously 
crucial as well. Unfortunately our local density can not be pinned down very well from N- 
body considerations alone since the simulations do not include baryons. The value of the 
local density used in most dark matter detection studies is p© ~ 0.3 GeV/cm^ [41,42], 
however this is based upon an older understanding of potential DM density profiles. A more 
recent study that has attempted to fix this value more carefully and construct a probability 
distribution for Pq, finds a canonical value of p© ~ 0.4 GeV/cm^ that can vary by up to a 
factor of ~ 10, but is likely to be less not more than a factor of 2 [43]. While the N-body 
simulations can not constrain this value very well, the recent study based on the Via Lactea 
II billion particle simulation finds a value consistent with this, albeit with a possible large 



(^''The Moore profile has been attributed to several values of a,/?, 7. The original profile, often referred to 
as M99, is very cuspy and has 7 = 1.5. Sometimes this is used in dark matter calculations and it artificially 
inflates the rate of annihilations in the center of the Galaxy. While there is no longer evidence for this cuspy 
of profile, a more modern study by Moore et al. [33] found a best fit with 7 = 1.2 which is usually denoted 
as Moore2004. 

(^'This is sometimes referred to as a Merritt profile in the literature [37]. 
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Figure 6: rp{r) for the different DM profiles. When appropriate, we also plot the regularized 
profile, with = 100 pc. Thick lines (from the bottom up): Einasto profile with a = 
0.2,0.17,0.12. Thin lines: isothermal, NFW, Moore (2004). The dashed lines represent the 
effect of the regularization described in the text. 

variation due to local clumping of dark matter in their simulation [44]. For our purposes 
we will use a value of Pq ~ 0.3 GeV/cm^ so our results can be easily compared with those 
of other groups. Nonetheless, but one should keep in mind this potential uncertainty when 
examining our results. 

Finally, the cuspy profiles which diverge as r — > 0, should be regularized. Similar to [45], 
we regularize the NFW [30] and Moore [32] profiles, by assuming at small r < Tc = 100 pc, 
the DM density takes the form 



p(r < Tc) , , . f TTr\ ^ . 2 f^f\ 
\ . ' = A1 + A2 smc — + A3 smc^ — 
p[rc) \rcj \rcj 



(3.4) 



The coefficients Ai are determined by demanding that p(r) and it's first derivative are con- 
tinuous and that the regularization does not alter the overall DM mass in the Galaxy. As 
we show later, the DM signal is rather sensitive to Vc- 

In Figure |6] we plot the various DM profiles as a function of the radial distance from the 
Galactic center (GC). When appropriate, we also show the regularized profile. 



3.2 Positrons 
3.2.1 Propagation 

To calculate the positron flux at the Earth one needs to understand how positrons propagate 
through our Galaxy. Due to our limited understanding of the latter, we wish to use a 
simplified model which employs a minimal set of assumptions. Perhaps the simplest model is 
the so called leaky box model which assumes a free homogeneous diffusion of charged particles 
within the Galactic disk. Since in its simplest form the model does not take particle cooling 
into account and since DM is not homogeneously distributed in the Galaxy, it is insufficient 
for calculating the DM signal. However, it is useful for calculating backgrounds and we will 



return to this model in Section 3.2.2 and the Appendix. 
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Below we adopt a widely studied diffusion model which does take cooling into account. 
Ignoring other effects in the propagation (such as convection and re-acceleration) is consistent 
for positrons at energies above ~ 10 GeV [28]. Moreover, the popularity of this model 
among the DM community allows one to compare our results with those of others, thereby 
emphasizing the new particle-physics features of our study. Where appropriate, we will 
remark on the uncertainties both in our theoretical understanding of the physical mechanisms 
involved and in the experimental tuning of the model [29]. 

To this end, the diffusion-loss equation in the steady state regime takes the form [46,47] 

- K{E)V^n,+ {E,x) - ^ {b{E)n,+ {E,x)) = Qe+{E,x). (3.5) 

Here ne+ {E, x) is the positron number density per unit energy which is related to the positron 
flux through $e+ = Ve+ne+ /Air. Qe+{E,x) is the source term for the positrons which for the 



background will be discussed in section 3.2.2 while for the DM source is given by 



Finally 

is the energy loss coefficient due to Inverse Compton Scattering (ICS) and Synchrotron Ra- 
diation, with te = 10^^ s. Here and below, E denotes the kinetic energy of the corresponding 
particle. The diffusive region is assumed to be a fiat cylinder parameterized by (r, z, 6), with 
z e [—L, L] and r e [0, i?], = 20 kpc. The positron number density is taken to vanish at 
the boundaries and particles may propagate freely outside of it. The solar system is located 
at (r0,Z0,6'0) = (8.5 kpc, 0,0) . 

Diffusion arises due to the interactions of charged particles with the galactic magnetic 
field inhomogeneities. Such interactions produce stable and unstable spallation products that 
may be used to extract the height of the diffusion region, L, and the diffusion coefficient 
K{E) which is usually taken to be of the form 

K{E) = Ko(3{n/GVY. (3.8) 

Here (3 = v/c and TZ = pc/eZ is the rigidity of the particle. For electrons and positrons, 
eTZ ^ E. The most stringent constraint on the above parameters comes from the B/C 
measurements of the HEAO-3 experiment [48] and other balloon experiments, the most 
recent being ATIC-2 [49]. The range for Kq, 6 and L has been studied e.g. in [27,50]. We 
adopt the parameters of [51], shown in Table [3] The (MIN, MAX) are taken to extremize 
the antiproton flux, while the (Ml, M2) extremize the positron flux. However, as it will 
be shown below, the differences between MIN and M2 are small for the positrons, while 
MIN encompasses a larger spread for the antiprotons. These sets of parameters are used 
as representatives of the uncertainties in the diffusion model. It is important however to 
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Model 


6 


Ko [kpcVMyr] 


L [kpc] 


Vc [km/s] 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX (Ml) 


0.46 


0.0765 


15 


5 


M2 


0.55 


0.00595 


1 


8 



Table 3: Parameters for the cosmic ray propagation model, compatible with measurements 
of B/C. The names correspond to minimum, medium or maximum primary antiproton and 
positron fluxes. The convective wind parameter Vc is only relevant to the propagation of 



antiprotons, as discussed in section 3.3 



note that these parameters are model dependent and their spread may be larger at high 
energies. For example, 6 is expected to change with energy, because an extrapolation of 
6 ~ 0.6 to higher energies predicts large anisotropics already at 10^^ eV, contradicting 
~[52]. 



observation; 



(3) 



Eq. (3.5) sets a natural energy-dependent diffusion length scale [29,53] 

V4i^o(tB 

where tE = te(1 



Xd{E, Es) = ^AKo{tE - tEs), (3.9) 
S)~^{E/GeVY~^ is the energy-dependent pseudo time. Using the models 



given in Table [3} we plot A d for detected electron with energy of 20 GeV, as a function of 
the injection energy in Fig. [7j A fundamental difference exists between the background and 
the DM contributions to the spectrum: while the injection spectrum of the former peaks at 
low energies (see eq. (3.14)), the one of the DM does not, as can be seen in Figs. 3]|5, It then 
follows that while background electrons travel short distances, of order Ad < 2 — 3 kpc, the 
bulk of the electrons originating from DM annihilations travel larger distances and may be 
created closer to the center of the Galaxy where the DM density is larger. We will return to 
this point shortly. 

A semi-analytic solution to the diffusion equation above is found to take the form [29] 



$e+(^,r0) 



ATrb{E) 



dE'Q,+ {E',rQ)I{XD{E,E')). 



where 



I{^d)= ^ Jo(Cn'^0/i?)sin(m7r/2)exp 



n,m=l 



\2l) 




Rr. 



and 



R 



rdr I dz JoiCn^ / R) sin {rmT{z + L)/ 2 L) 



P@ 



(3.10) 



(3.11) 



(3.12) 



'^'^•'it is possible that such a change is indicated in the ATIC-2 measurements [49] (which are not taken 
into account in [50] ) , although the large statistical uncertainties and possible normalization problems do not 
allow for a decisive conclusion. 
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Figure 7: Left: The effective diffusion length, X^, for the propagations models we study 
as a function of the injection spectrum for a 20 GeV electron detected at the Earth. Right: 
/(Ad) for L = 1, 4, 15 kpc and the Einasto profile with a = (0.12, 0.17, 0.20). 



Here Ji are the i-th order Bessel functions of the first kind and (n is the n-th zero of Jq. 
We plot /(Ad) for the Einasto profiles with a = (0.12,0.17,0.20) and in Figure [T] for the 
three different values of L considered here. As can be seen in the plot, /(Ad) may provide an 
enhancement to the flux if Ad is sufficiently large. Since Ad is largest for positrons arriving to 
the Earth with low energy, the spectrum tends to be softer in those cases where the positrons 
propagate for significant distances. This, in turn, weakens any feature that may exists at 
the high end of the spectrum. As suggested in Fig. [7], the effect is most significant for the 
MAX model. Consequently, as we shall see in Sect. |4], the MAX model fits less well to the 
ATIC data, but on the other hand requires a smaller cross-section to fit the PAMELA data, 
therefore predicting less photons from the GC. Conversely, the MIN model fits better to the 
ATIC data but predicts a larger cross-section and therefore a larger number of photons. This 
tension, in part, will be responsible for excluding significant parts of the parameter space in 
these DM annihilation models that try to fit PAMELA and ATIC simultaneously. 



3.2.2 Experiments 

The three relevant experiments for the positron study are PAMELA [1], ATIC-2 [3] and 
PPB-BETS [2]. The former recently provided the positron-to-electron flux ratio up to 100 
GeV while the latter two experiments measure the positron plus electron flux up to energies 
of a few TeV. 

For the purpose of testing the predictions of the theory, we need an estimation of the 
backgrounds. To this end, we will attempt to be as data-driven as possible. The two relevant 
backgrounds needed are the positron fiux in the energy range between 10 to 100 GeV and the 
e~^+e~ fiux at energies up to 1 TeV. Starting with the latter, the spectrum has been measured 
rather accurately. In the absence of the ATIC/PPB-BETS bump, which is apparent above 
~ 100 GeV, the data fits a power-law spectrum 

^e++eAE) = A, (y^Y cm-^ s-1 sr-1 GeV-^ (3.13) 
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For finding a best fit to Ai and 71 we used the measurements of ATIC-2 and HEAT [54] in the 
energy range 20 — 100 GeV. The lower cutoff ensures that uncertainties from solar modulation 
are insignificant while the upper prevents probing the bump. We find Ai = 0.046 ± 0.015 
and 7i = —3.30 ± 0.08. This result is in agreement with the preliminary analysis presented 
by the PAMELA collaboration [55]. Taking the best fit at face value for our background 
is incorrect because it is impossible to disentangle the signal from background even at low 
energies. Therefore, when fitting together with the DM signal, we allow both Ai and 71 to 
float around their central values given above. 

No precise measurements of the positron flux exist at energies above 10 GeV. To calcu- 
late the background, we therefore use the leaky box approximation which is well suited for 
energies < 100 GeV and above > 10 GeV^^^ We derive the prediction in Appendix [a| and 



simply quote the result here 

/ T? \ -2.84-5±0.02 

^e+{E) ~ 4 X 10"^ ( j cm-2 s~^ sr-i GeV"^ (3.14) 

which is in agreement with [28]. The parameter 5 is typically of order 5 ~ 0.6 but for 
the propagation models considered above, we vary it according to the values of Table [3j 
As described in the Appendix, the normalization of the flux suffers from large uncertainties 
which arise from the assumptions in the theory and uncertainties in the measurements of the 
nuclear spallation cross-sections. Thus when fitting the data together with the DM signal, 
we allow the normalization to float around the above value. The spectral index is instead 
kept fixed. 

3.3 Antiprotons 
3.3.1 Propagation 

We describe the propagation of p at equilibrium with a diffusion equation similar to that of 
the positrons [56] 

- K{E)V^np{E, x) + -^ {Vc{z)np{E, x)) + 2h6{z)T,^nnp{E, x) = Qp{E, x) (3.15) 
where np{T, r) is the number density of antiprotons per unit energy. As in the case of 



positrons, K(E) is the diffusion coefficient given in Eq. (3.8). Energy loss, however, is 



negligible for antiprotons due to their mass and we therefore do not consider it. The last 



two terms on the l.h.s. of eq. (3.15) correspond to convective wind and interaction with 
the Interstellar Medium (ISM) in the galactic plane respectively. The convective wind is 
assumed to be of the form Vc{z) = sign{z)Vc and is directed outwards from the galactic 
plane. The annihilation width between antiprotons and protons is given by 

Tann = {uu + 4'/=^nHe)a|°%, (3.16) 



'^^•'We thank Eli Waxman and Boaz Katz for drawing our attention to this point. 
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where nH,He — 1 cm ^ are the hydrogen and hehum number density and Vp is the velocity of 
the antiprotons. A parameterization for be found in [56,57]. Both the annihilation 

and convective terms can be neglected at energies above > 10 GeV. Nevertheless we keep 
these terms for completeness. 

As in the positron case, propagation takes place within the disk of half-height L and 
radius i? = 20 kpc. In order to be able to fit the data both for positrons and antiprotons 
and since we would like to quantify the uncertainties, we use the same models described in 
Table [3] An semi- analytical solution to the flux $p = VpUp/An is given by [58] 



ME, re 



Ait 



Qp R{E) 



(3.17) 



where as before. 



Qp 



dE' 



2 mJ^ 

k 



(3.18) 



R{E) encodes the astrophysical propagation information and was found to take the form [51], 



R{E) = Y,MCnrQ/R)ew 



n=l 



VcL 



2K{E) 



yn{L) 



yl„sinh(5„L/2) 



(3.19) 



with, 

yn{L) 



(3.20) 



X 



R 



rdr / dzJo{(nr/R) sinh(5'n(L — z)/2)exp 
Jo 

2/ira„„ + Vc + K{E)Sn coth(S„L/2), 



Vc{L - z) 
2K{E) 



p{r, z) 
Pq 



Sr, = ^{Vc/K{E)y + 4aR^. 



In Figure [8] we plot this function for the the Einasto profiles with a = (0.12,0.17,0.20) and 
for the propagation models of Table [3] 



3.3.2 Experiments 

The relevant experimental data is the recent PAMELA measurement for the antiproton-to- 
proton flux ratio [8] up to the energy of lOOGeV. 

The proton flux, has been well measured by the BESS [59,60] and AMS-01 [61] experi- 
ments and the best fit to the data (for demodulated protons) is of the form [62]: 

= AP^'^n-'"' cm-2 s-i sr-i (GeV/n)-' (3.21) 
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Figure 8: R{E) for the propagation models we consider and the Einasto profile with a = 
(0.12,0.17,0.20). 

where as before TZ is the rigidity and (3 is the velocity. We use the following values for the 
coefficients 



n 


A 


Pi 


P2 


< 20 GV 


1.94 ±0.13 


0.7 ±0.52 


2.76 ±0.03 


> 20 GV 


2.42 ±0.18 





2.84 ±0.02 



in agreement with [62]. 

The antiproton fiux is known to come primarily from cosmic ray protons interacting with 
the ISM. Unfortunately the fiux generated by this process does not have a characteristic 
power law shape (as many other astrophysical processes) in the range of energy that we ex- 
amine. To calculate the approximate shape for the background we need to take into account 
the spallation cross section given in [62], and then propagate the products. Fortunately, 
as was found in [62], changing the model of propagation does not significantly change the 
shape of the antiproton fiux but does change the normalization. We therefore use the results 
of [62], parametrized as follows [26], 

log^g $p = -1.356 ± 0.114r - 0.645r2 - 0.614r3 + 0.02r^ ± 0.168r^ - 0.038r^ (3.23) 

where r = Xog^Q^E / GeV) . To fit to the PAMELA data one needs to take into account 
solar modulation [63] . We modulate the background according to the PAMELA modulation 
potential = 0.58 GV and allow the normalization to fioat. 



3.4 Photons 
3.4.1 Propagation 

Photons from dark matter are one of the cleanest channels possible to study from the point 
of view of propagation. This is because once produced, photons freely travel to Earth and 
thus are insensitive to the propagation parameters that complicate studies of e"*" and p. Since 
photons are not affected by propagation, the fiux at the Earth depends only on the injection 
spectrum and on how many annihilations occur. While this simplifies the calculation of 
the photon fiux, it also introduces a large amount of uncertainty. Unlike positrons and 
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antiprotons which reach us typically from distances less than r©, photons from the galactic 
center, where the dark matter density is the highest, can always reach us. N-body simulations 



can not resolve distances less than approximately 100-200 pc (as discussed in 3.1), and thus 
there is a fair amount of uncertainty in the DM density profile at the galactic center. 



The differential photon flux (in units of cm s sr GeV ) is given by. 



Air 



Q., J An 



(3.24) 



As before, is the photon source encoding the particle physics contribution to the flux. 



lp(r) 



2 m^2 

^ k 



Em, 



'~dE' 



(3.25) 



Much like I{E) and R{E) in the positron and anti-proton case, J encapsulates the astro- 
physics dependence and is given by. 



J 



An 



dn 



dl f p{r) 



line-of-sight '"O \ P& 



(3.26) 



where An is the opening solid angle for a given experiment. The flux received in a given 
direction is proportional to the integral over the line-of-sight of the DM density squared. By 
writing the photon flux in this factorized form we simply need to compute the J for a given 
experiment and profile. 

In principle, for cuspy profiles the integral in Eq. 3.26 could diverge. For the profiles 
that we consider J in fact does not diverge. However, as discussed in 3J^, depending on the 
assumed regulation of the DM density inside of C(IOO) pc, the actual value of J computed 
can vary significantly as a function of in Eq. 3.4 We will quantify this uncertainty further 



in Section 4.2| when we show the dependence of photon flux on the choice of the dark matter 
density profile. 



3.4.2 Experiments 

The experiments that we will be primarily concerned with are those that measure high energy 
gamma rays. Specifically we will focus on the HESS experiment which is the most sensitive 
to high mass DM annihilation. The FERMI-LAT experiment will also be of importance in 
the near future and we discuss it further in Section [5] If one is looking for signals of DM 
annihilation into photons, the best place to look is in regions of potential high DM density 
with small astrophysical backgrounds. For instance high mass/light ratio dwarf spheroidal 
galaxies [64] or subhalos of our galaxy [65] are both prime candidates for a clean signal 
of DM annihilations. However, for the purpose of constraining models of DM, it is more 
important to observe regions that have effectively a large J for sufficiently long time. For 
the HESS experiment this is found in two regions, the Galactic Center (GC) [66] and the 
Galactic Ridge (GR) [67]. 
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The GC as studied by HESS has a AQ = 10 ^. To calculate J for the Galactic Center 
one can write 

Jdn = J^^^^27rd{cosi;') , (3.27) 
An = 27r(l - cos ^Z-) , r = ^^r% + P - 2lrQ cosip'. (3.28) 

The GC was observed by HESS for 48.7 hours during 2004, and is cataloged as J1745-290. 

The GR is a larger region that includes the GC. It is defined in galactic latitude and 
longitude as the region \l\ < 0.8° and |6| < 0.3°. The GR has several known suspected 
point sources including J1745-290 that the HESS experiment modeled in order to subtract 
it in [67]. The experiment fits a point source to two specific regions, J1745-290 and another 
source, GO. 9+1, that lie within the Ridge. The best fit for these sources was then subtracted 
and the differential flux dN/dE was given for the remaining portion which defines the Ridge. 
Unfortunately the dN/dE for each galactic latitude and longitude has not been published 
and therefore one can not mimic their subtraction scheme. However, since the hypothesized 
point sources are the dominant sources of photons we can excise these entirely without 
introducing a large amount of uncertainty into the calculation of J . Additionally, beyond 
these partial subtractions, a further subtraction to reduce systematics is also performed. 
The HESS experiment observes the regions 0.8° < |6| < 1.5° and uniformly subtracts this 
"background" from the GR data. We will perform this same subtraction when plotting the 
observed photon flux. To calculate the value of J for the GR (which is a rectangle not 
a circle) we use a slightly different form of dVt. For a region defined by longitudinal and 
latitudinal boundaries \9\ < A/, < A6, we have, 

jdn = A jf' dO' /o""^' d{sm^') , (3.29) 
An = AAlAh , r = ^Jr% + P - 2lrQ cos ^j' cos 6'. (3.30) 

To finish our discussion of the experimental inputs for photons we will tabulate the J 
calculated for the various experiments and profiles in Table |4j 



3.5 Neutrinos 
3.5.1 Propagation 

Dark Matter annihilations can also produce neutrinos. After being produced, neutrinos 
propagate till the Earth where they are detected. The standard strategy to detect these 
neutrinos is to look at those that convert in rock nearby the detectors and observe the 
charged leptons, in particular muons. Therefore the quantity actually measured is the muon 
flux, that can be related to the neutrino flux at production by 

d^^ /"""^P p da,N{E^/E,) 

— = R{E,) dE, Y: ^P^-^ E ^^-'^-) (3-31) 
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Location 


Profile 


J 


Galactic Center 
{An = 2 ■ 10-5) 


Cored isothermal 
NFW 
Moore 
Einasto a = 0.17 
Einasto a = 0.12 
Einasto a = 0.2 


13.6 
4076 - 10170 
13128 - 51388 

6610 

65500 

2306 


Galactic Ridge 


Cored isothermal 
NFW 
Moore 
Einasto a = 0.17 
Einasto a = 0.12 
Einasto a = 0.2 


0.02 

1295 - 1541 
3836 - 5653 
1614 
10886 
602 



Table 4: The values for J for the various galactic regions we consider. For the NFW and 
Moore profiles in the Milky way, we show the range for J for < Tc < 0.1 kpc. 



where d^^^JdE is the flavor speciflc neutrino flux, the sum is taken over all the leptonic 
flavors and the term Pj^^ comes from the fact that neutrinos oscillate during their trave 



(5) 



The muon production cross section in neutrino nucleon interactions is denoted by a^^ [68,69], 
while pp^n are the proton and neutron number densities (taken to be equal for standard rock). 
Finally R{E) is the muon range [42], i.e. the distance traveled by a muon of energy E before 
stopping or before getting below the detection energy threshold of the experiment. 



3.5.2 Experimental Status 

Good places to look for neutrinos produced in DM annihilations are the GC (like in the case 
of photons) but also the Sun and the Earth, where DM may be trapped in the gravitational 
field and their density may grow large enough to allow for a sizable annihilation rate. Another 
signal that one can look for is the total diffuse neutrino flux [70]. 

The flux coming from the Sun and the Earth has been recently re-investigated in light 
of the PAMELA and ATIC [71]. In particular the authors of [72] have found that for DM 
annihilating mainly into leptonic flnal states, the neutrino flux from the Sun and the Earth 
is out of reach of IceCube unless the DM annihilates directly into a pair of neutrinos. In our 
case IceCube can still be of relevance if the mass is high enough in the multi-GeV range 
where hadronic decays are also allowed, but in this case the bounds on the photon flux will 
be also very strong as it will be shown in Sect. |4j For this reason here we will only discuss 
the neutrino flux coming from DM annihilations in the Galactic Center. 

(^^Pi^^ is determined only by the mixing angles since the baseline is very long. For simplicity we have 
taken 6*13 — 0. 
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In this case the present bounds comes mostly from detectors located in the northern 
hemisphere, since there the Galactic Center is below the horizon most of the time. In 
particular the current best limit is from SuperKamiokande [73], which looks at the up-going 
muons produced in the rocks below the detector. The collaboration reports an upper bound 
on the total muon flux above a threshold of 1.6 GeV, as a function of the half-opening angle 
of a cone pointing towards the center of the Galaxy. As we show in the next Section, this 
bound is not powerful enough to give any appreciable constraints on the parameters of our 
module. 

In the future this limit will be improved by Antares, a neutrino telescope in deep Mediter- 
ranean waters, which started taking data in 2007, and on a longer timescale by Km3Net and 
Megaton-size neutrino detectors like Hyper-Kamiokande. 



4 Results 

In this section we combine the particle physics inputs that were calculated in Section [2] and 
the propagation methods that were discussed in Section [3] to calculate the fluxes for various 
experiments. We are primarily interested in answering the following questions: 

1. Given our particle physics framework, what is the preferred parameter space that can 
explain PAMELA and ATIC/PPB-BETS? 

2. Given this parameter space, what are the regions that are compatible with existing 
searches for 7's and z/'s? 

To answer these questions we need to define what we mean by the preferred parameter 
space and how we bound it. We implement a function in order to fit a given data set 
with our particle physics and background parameters (the relevant set is shown in Table [s]). 
To understand the allowed regions and how the parameters affect the fluxes we perform 
the following procedure. We pick a parameter of interest and we marginalize over all the 
other parameters. We then extract the fluxes as a function of the chosen parameter. In this 
section we will only marginalize a function for the PAMELA and ATIC/PPB-BETS data. 
We then plot the resulting fluxes as a function of this parameter both for PAMELA, ATIC 
and PPB-BETS, and for HESS and SuperK data, where applicable. The allowed parameter 
regions are those points with good flts that do not conflict with the HESS and SuperK data. 

Below, we separate the parameters that we are interested in into two categories: particle 



physics discussed in Section |4.1[ and astrophysics discussed in Section |4.2[ By examining 
these parameters independently, we study their influence on the fluxes, and the uncertainties 
in our predictions. We postpone the flts to all data and their implications to Section |5| 

4.1 Particle Physics 

4.1.1 Light Gauge Boson Mass 

In this section we isolate the effects of varying the light gauge boson mass, m^, on the fluxes. 
To understand how the particle physics parameter influences the results, we flx both the 
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Background Parameters 


Ne+, Ne++e- , Np , 7e++e-, ,7^ 


Particle Physics Parameters 


m^, , {av) , i?sM, odm 



Table 5: A summary of the parameters used for fitting the data. 



(GeV) 


(TeV) 


10^^ X (cTf) (cm^s ^) 


(iV,7)e++e- 




xVdof 


0.2 


0.96 


0.94 


(0.6,-3.21) 


0.7 


1.9 


0.5 


1.54 


2.70 


(0.6,-3.20) 


0.6 


2.1 


1.2 


1.46 


2.51 


(0.6,-3.19) 


0.6 


2.1 


8 


2.3 


5.41 


(0.3,-3.06) 


0.2 


2.7 


15 


3.5 


9.17 


(0.3,-3.06) 


0.3 


2.9 



Table 6: The best fit values for the plots in Fig. |9j 7e++e- is the best fit value for the spectral 
index of the background electron plus positron fiux. A^^ is the fraction of the normalization 
found for the best fit background without a DM signal. These normalizations are found in 
Section |3| Finally, the SM ratios, Rsm are all smaller than 1%. and so are not shown. 

dark matter profile and propagation model and marginalize over the rest of the parameters. 
We choose the Einasto profile with a = 0.17 and the MED propagation model. 

In Fig. |9| we plot the best fits for PAMELA, ATIC/PPB-BETS as well as the bounds 
from neutrinos and HESS experimental data for the GC and OR in Figure [9j In Table |6] 
we give the best fit values for these plots together with the x^/dof. The number of dof's 
that go into our fit is 33, from assuming 40 independent measurement points (coming from 
PAMELA, ATIC and PPB-BETS) and 7 variables. As one can see, the x^/dof is typically 
of order 2 which for 33 dof's is a rather poor fit. We trace this fact to the two features in 
the ATIC data: The plateau above 100 GeV and a bump at around 500 GeV which cannot 
be explained together with only one DM state. It remains to be seen if these two features 
in the spectrum survive future measurements. In the mean time, we view the values of 
as no more than an indicator for comparing different models. 

The values plotted for are discussed in Section |2| From Table |6] we see that small m,^ 
fit the data better. This is because as increases the positron spectrum becomes softer due 
to new kinematically accessible decay channels. As increases dN/dx increasingly deviates 
from the relatively fiat spectrum when only annihilates into e"*" + e~. Compared to the 
fiat dN/dx, to fit PAMELA and ATIC/PPB-BETS for higher m^, one needs a combination 
of boost factor and increase in mass to attempt to make up for the softer spectrum. 

For the case of 8 and 15 GeV, the large DM mass (and therefore the poor fits) is driven 
by the more constraining antiproton fiux. Indeed, for above twice the proton mass, the 
PAMELA antiproton-to-proton spectrum becomes an even stronger constraint. In particular, 
for = 8 and 15 GeV, the DM mass must be sufficiently large in order to avoid generating 
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PAMELA Poxiirons 



PAMELA Aiiiiprotoiis 




Figure 9: Best fit to the PAMELA, ATIC and PPB-BETS experiments for Einasto profile 
with a = 0.17, propagation model MED fixed and = 0.2, 0.5, 1.2, 8, 15 GeV. The predic- 
tions for neutrinos and photons are shown in the last three plots overlayed with the SuperK 
and HESS measurements respectively. Around each bin of the HESS data, we indicate the 
Icr (solid) error bar. For the GR we also indicate the 3cr (dasshed) error bar. The best fit 
parameters are shown in Table [6j 
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too many antiprotons at energies within the PAMELA reach. This is indeed apparent in 
Fig. [9] which shows a large bump in the antiproton spectrum just above the PAMELA 60 
GeV bin. Whenever the antiproton constraint pushes the mass too high, denying a good fit 
to the positron data, the overall fit for these masses is significantly worse. 

It is interesting to note that a 2 — 3 TeV DM particles is sufficient in order to avoid the 
antiproton bound in the case of 8 GeV light gauge field. This is in contrast to a 10 TeV DM 
required in the case of direct W and Z decays as in the SM [74]. A simple way to understand 
this, is to notice that the anti-proton spectrum is (almost) identical in both cases when 
viewed in the gauge boson restframe. The spectrum must then be boosted to the lab frame, 
which introduces a boost factor of order m^/m^. In both the SM and with light gauge fields, 
this factor is of order 10^ which is the minimum required in order to avoid the PAMELA 
antiproton bound, effectively increasing the QCD scale above 100 GeV. Thus we conclude 
that with light gauge fields the bound is easy to overcome with, m^/TeV ~ m^/(10 GeV). 

The photon predictions for the HESS measurements from the GC and GR are shown in 
the last two plots. It is apparent that the GR is more constraining than the GC. Moreover, 
we do not include photons from ICS which may strengthen the bound. Heavy vector bosons, 
at or above 1.2 GeV are excluded by HESS. On the other hand, the 200 and 500 MeV vector 
bosons are marginally within 3a error bars of the measured data. Clearly, once a background 
is added to the fit the data, the tension becomes larger, as the background is required to 
be sufficiently large to allow for a good fit to the data points above where a DM signal 
is absent. We therefore conclude that some further suppression (e.g. local boost factor or 
shallower DM profile) is required to render these models viable. We return to this issue and 
quantify the suppression in Section |5j 

Finally we note that the Neutrino bound is sufficiently weak to evade and therefore adds 
no further constrain on the model. 



4.1.2 Hidden Sector Shower 

Let us now study the effects of showering in the dark sector. As in the case where we 
isolated the effects of m^, we fix an Einasto profile with a = 0.17 and MED propagation, 
then marginalize over the other parameters and plot the fluxes for different values of aoAf- 
Speciflcally to illustrate the effects of the shower we flnd the best flt for a^M = 10~^ and 
then keeping the same parameters but varying the gauge coupling we plot the resulting 
fluxes in FigjlOj To illustrate the effect of showering on the antiprotons, we repeat the same 
procedure with = 8 GeV, showing only the antiproton flux ratio. 

There are several effects worth noting. The flrst, is that showering in the DM sector 
implies a softer electron, positron and antiproton spectra. As discussed in the previous 
section, a softer spectrum flts ATIC less well. Moreover, if antiprotons are produced the 
softening of their spectrum forces to be higher. On the other hand, softening of the 
spectrum implies more electrons and positrons at the PAMELA energies, which explains the 



increase in flux for larger odm as can be seen in the flrst plot of Fig. [TO 

On the other hand, as we will see in Section |5j models where showering in the DM sector 
is signiflcant, tend to flt the data better when taking the HESS measurements into account. 
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Figure 10: Effects of dark sector sfiowering. Plotted are tlie PAMELA, ATIC/PPB-BETS 
and HESS GR data for models witli Einasto profile witfi a = 0.17, MED propagation model, 
= 200 MeV, = 1000 GeV. and odm = 0.001, 0.01, 0.04, 0.1. These parameters are the 
best fit for the qdm = 0.001 case. In this model antiprotons are not produced. Instead we 
show the effect on antiprotons for the case of = 8 GeV. Around each bin of the HESS 
data, we indicate the la (solid) and Sex (dasshed) error bar. 



The reason for this is that the photon spectrum is softened as well due to the showering. This 
is apparent in the last plot of Figure 10 where we see that larger a^u implies less photons. 
Since the strongest constraints come from photons, it is preferred (for minimizing the x^) to 
allow for softening of both the photons and electrons, so that the ATIC fit is poorer but the 
photon fit is better. 



4.1.3 SM Ratio and Direct Decays 

We now study the implications of DM directly annihilating into SM particles. A similar 
study was presented in [15,74]. Here we confirm part of their results, while emphasizing the 
difference between the predictions of direct annihilations versus those through light vector 
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fields. 

To this end, we concentrate on two possibilities, namely direct couplings to SM gauge 
bosons, and direct couplings to e^. From the point of view of the PAMELA results, the latter 
is preferable as it does not incorporate any hardronic decays and therefore no anomalous 
antiproton flux. On the other hand the former is guaranteed to exist at some level if the (p 
can mix with the SM vector bosons. 

To emphasize the difference between 2^2 and 2^4 decays, we compare the predic- 
tions of direct anihilations to the case of = 200 MeV, for which only production is 



kinematically accessible. As was discussed in subsection 4.1.1, this is also the model which 
best fits the data if shower is not included. The best fit parameters are shown in the first 
line of Table ^ The value for the best fit of the direct decay is found to be x/dof = 2.6 



which is worse than that of the 2 — 4 case. The best fit parameters are. 



m 



X 



680 GeV, {av) = 6.34 x lO^^s ^^^2^ j^^^ ^ 2%, (4.1) 
(iV, 7),++e- = (0.9, -3.27), = 1. (4.2) 



As in previous cases, Ni denotes the fraction of the normalization found for the backgrounds, 
while 7e++e- is the spectral index for the correspond background flux. 



In Fig 11 we plot the two models together with the PAMELA positrons, ATIC/PPB- 
BETS and HESS measurements. While the fit to the ATIC/PPB-BETS data is better, it is 
clear that direct coupling to does not fit the PAMELA data well. The reason for this can 
be traced back to the injection spectrum for the line shown in Fig. |4j Since the decay is 
2 — >■ 2, the ATIC data dictates a low DM mass around the bump. For such a low mass, the 
hard form of the spectrum is then insufficient to explain the PAMELA anomaly. 

To partly compensate for the injection spectrum, a large cross-section, is required. This, 
together with the very hard spectrum explains the large number of photons predicted in such 
a case, as shown in Fig. [TTj Clearly, this model is in contradiction with the HESS data and 
is therefore excluded [15]. 

Next we would like to understand the extent to which the DM particle can couple to 
the SM gauge vector bosons. A known constraint is that the massive SM gauge fields decay 
into hadronic states and may therefore produce an unacceptable excess of antiprotons. It is 



interesting to quantify this statement. We do that in Figure 12 where we plot 2D confidence 
level contours for the allowed regions of parameter space with = 200 MeV and Einasto 
a = 0.2 profile, as a function of and Rsu- We see that while the best fit value prefers 
5% of SM annihilations, at la one can have as much as 25%. While in this plot we did 
not attempt to fit the HESS data, we have checked that models which fit the HESS data 
without difficulty, such as Einasto profile with a = 0.2, do not change this result significantly 
(see Sect. |5]). Conversely, models that are only marginally consistent with HESS, such as 
the Einasto profile with a = 0.17, reduces the parameter space almost entirely. For model 
building this result implies that the DM particle must couple only weakly to the SM. This can 
be implemented either by coupling x to SM singlets, or otherwise making odm sufficiently 
large. 
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Figure 11: A comparison between an e^e~ annihilation channel and the best fit result for 
annihilation through two 200 MeV light state which each decays into e^e~. The details of 
the 2 — s> 4 model is shown in Table [6j The fit done with the Einasto profile with a = 0.17 
and propagation model MED fixed. The antiproton and neutrino spectrum is not shown 
since those are not produced due to kinematics. The predictions for photons are shown in 
the last two plots overlayed with the HESS measurements. Around each bin of the HESS 
data, we indicate the la (solid) and 3cr (dasshed) error bar. The best fit parameters are 
shown in Table |6l 

4.2 Astrophysics Uncertainties 
4.2.1 Profile 

We now study the dependence of the predictions on the DM profile. As discussed in Sec- 
tion 3.1[ current N-body simulations do not allow us to pin-point the precise DM profile. 



The main difficulty for these simulations is the resolution, which does not allow one to probe 
the DM distribution within ~ 100 pc from the GC. Moreover, baryons are not incorporated 
in simulations, while they may play important roles in the GC. 

The effective diffusion scale for electrons is smaller than the distance of the solar system 
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Figure 12: A confidence level contour plot for m,^ = 200, and Einasto a = 0.2 as a function 
of and Rsu- The contours correspond to the 68%, 95% and 99% confidence levels. 



from the GC and therefore their fiux is not sensitive to the large uncertainties in the inner DM 
profile. On the other hand, photons do not diffuse and therefore most of them come from the 
center where the bulk of the DM lies. This then allows one to probe and constrain theories 
of DM in conjunction with the DM profiles that are extracted from N-body simulations. 

To study the profile dependence we take = 500 MeV, and scan over the six profiles: 
NFW [30], Moore [32], Isothermal Core [34] and Einasto [36] with a = (0.12,0.17,0.20), 
fitting the rest of the parameters to PAMELA and ATIC/PPB-BETS. Because we do not 
attempt to fit to the HESS data here, the fits are almost identical in their resulting param- 
eters. We find (on average). 



= 1540 GeV, 



{av) 
(0.6,-3.21) 



2.88 X 10"23 cm2, 



0.4 



-RsM 
■ 1. 



0, 



(4.3) 
(4.4) 



In Figure 13, we show the usual plots together with the predictions for the photons. The 
bands in these plots show the sensitivity of the NFW and Moore profiles to the distance 
< < 100 pc from the center of the Galaxy, below which we regularize the profiles as 



discussed in Section |3.1[ No such regularization is performed for the Einasto and Isothermal 
profiles since they do no diverge at the GC. 

The HESS plots in the figure demonstrate the strength of the constraint arising from 
the HESS data, the strongest coming from the GR. These essentially constrain J defined in 



Eq. (3.26). We stress that, as before, the expected photon spectrum shown does not include 



any background. Adding such a background is crucial in order to fit the measured data at 
energies above that of the DM. In such a case, the constraints are much stronger. We return 
to this issue in the next section. 
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Even without background we learn that both the Moore and Einasto with a = 0.12 are 
excluded by more than 3a while NEW and Einasto with a = 0.17 are above the data, but 
within the 3a error bars. As we shall see below, once background is added, these profiles are 
excluded unless some kind of suppression is in place, like the hidden sector shower, a local 
boost factor, etc.. 

It is important to note that independently of J, these photon predictions provide an 
conservative estimation of the signal for two reasons: (i). We have not taken ICS into 
account. Such radiation is certainly important at low energy, but could also be important 



at intermediate scales, (ii). As discussed in Section 3.4, for the GR, we do not include any 
photons coming from the GC. On the other hand, the GR measurements count photons from 
the center, after removing a dominat background source [67]. This is approximately a 10% 
effect. 

The annihilating DM scenario studied here is somewhat complimentary to the studies 
of [15]. We therefore conclude that the HESS measurements together with PAMELA and 
ATIC/PPB-BETS, strongly constrain the possible DM profiles in the case of annihilating DM 
scenarios. While not excluded, future experiments may strengthen the bounds considerably. 

4.2.2 Propagation 

Finally, let us consider the uncertainties arising from the propagation parameters. Here we 
stress again that our choice of propagation model already entails some uncertainties as our 
understanding of Cosmic ray propagation is not complete. Nevertheless, we show how our 



predictions vary as we change the parameters as discussed in Section 3.2.1 which were shown 
to span the possible positron and antiproton spectra. Such uncertainties have been studied 
extensively [28, 29, 75-77] so we only concentrate on our predictions. 

As in previous sections, we take our benchmark point, = 200 MeV and the Einasto 
a = 0.17 profile, and fit the rest of the parameters for the four propagation models, 
MIN/MED/MAX(M1)/M2. The best fit values of the parameters are shown in Table 



and the predictions are shown in Figure 14 Even though antiprotons are not produced 



from DM annihilations through the light gauge fields, some of the best fit models require a 
non-negligible SM couplings and consequently antiprotons and produced as can be seen in 
the plots. 



The tension anticipated in Section 3.2.1 is now apparent: because of the enhancement of 
the flux through I{Xd) for large Ad in the MED and MAX cases, the electron spectrum in 
that case is softer (recall that the diffusion length, A^), is larger for lower detected energies). 
Therefore, for these propagation models the feature in the high energy spectrum is less 
pronounced and so the fit for ATIC/PPB-BETS is not as good. On the other hand, due 
to the soft spectrum, a smaller cross-section is needed to fit PAMELA and therefore less 



photons are predicted for these models. This is consistent with Figure [I4| Moreover, in 
the MED/MAX models the spectral index for the background + e~ flux is harder to 
compensate for the bad fit to the ATIC bump. 

We find it is easier to evade the HESS constraints if propagation of positrons is closer to 
the MAX model described above, and in particular if the dependence of the escape time of 
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Figure 13: Best fit to the PAMELA, ATIC and PPB-BETS experiments for different DM 
profiles and with = 500 MeV and = 1.54 TeV. Propagation model is taken to be 
MED. No antiproton bound exists due to kinematics. The bands in the photon predictions 
come from regulating the NEW and Moore profiles with < < 100 pc. Around each 
bin of the HESS data, we indicate the la (solid) and 3a (dasshed) error bar. The best fit 
parameters are shown in Eqs. (4.3), (4.4). 
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Model 


(TeV) 


10^^ X ((Tf)(cm^s ^) 


(iV,7),++,- 




-RsM 


xVdof 


MIN 


0.90 


1.15 


(0.9,-3.28) 


1.9 


6% 


1.8 


MED 


0.96 


0.94 


(0.6,-3.21) 


0.7 


1% 


1.9 


MAX(Ml) 


1.38 


0.38 


(0.4,-3.10) 


1.4 


1% 


2.3 


M2 


2.3 


1.44 


(1,-3.33) 


0.9 


11% 


2.0 



Table 7: The best fit values different propagation models with m,^ = 200 MeV, shown in the 
plots of Fig. [14| 7e++e- is the best fit value for the spectral index of the background electron 



plus positron flux. Ni is the fraction of the normalization found for the best fit background 
without a DM signal. These normalizations are found in Section |3j 



Cosmic rays on the energy is weaker (smaller 6). Interestingly, it is likely that the behavior of 
the escaping time flattens somewhere below the knee at 10^^ eV [52]. Moreover, uncertainties 
in the spallation cross-sections may point towards a smaller 6. 



5 Implications and Future Directions 

In this section we discuss the implications for models that attempt to explain the PAMELA 
and ATIC/PPB-BETS excesses, based on the results shown in Section |4j In Section |4] we 
have shown by varying the particle physics and astrophysics parameters that only certain 
regions of parameter space can satisfy all the various experimental constraints. Specifically, 
the most difficult constraints arose from the HESS's measurement of the GC and GR regions. 
Indeed it was initially believed that hadronic activity from DM was needed to be suppressed 
in order to avoid creating an excess in the antiproton flux. Conversely, could be pushed 
to scales of order 10 TeV, to avoid the bound. As we have demonstrated in Section [4], one 
can tolerate antiprotons without having to raise significantly above the mass scale of the 
purported ATIC/PPB-BETS excesses. On the other hand, the real hadronic danger comes 
from 7r° decays which produce a significant amount of photons. By scanning over we 
have shown that one does in fact need a model that goes ultimately almost exclusively into 
purely leptonic final states. 

In several recent papers [15, 16] it was also argued, or indirectly demonstrated, that the 
high energy photons rule out models that try to explain both the PAMELA and ATIC/PPB- 
BETS excesses. This is not a conclusion that we agree with. These papers primarily in- 



vestigated the case where DM annihilated directly into leptons. As shown in Section 2A 
the dN/dx calculated in these cases is actually quite different, and a simple rescaling of 
the overall rate does not interpolate between the case of directly annihilating into leptons 
and annihilating through 0. To demonstrate that there exist viable models that satisfy the 
high energy constraints we show a best fit point from the scan over our parameters. In 



Figure 15 this point is shown where = 1 TeV, = 200 MeV, {av) = 8.2 ■ 10 cm /s 



and a^M = 0.04, for an Einasto profile with a = 0.2 and the MED propagation parame- 
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Figure 14: Plots for different propagation models, with = 200 MeV. The best fit param- 
eters used are given in Table [7j The neutrino spectrum is not shown since non are produced 
due to kinematics. The predictions for photons are shown in the last two plots overlayed 
with the HESS measurements. Around each bin of the HESS data, we indicate the la (solid) 
and 3(7 (dasshed) error bar. 
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Figure 15: The best fit including HESS/PAMELA/ATIC/PPB-BETS is plotted in blue(lower 
blue curve is signal only in the HESS plots), which is for = 1 TeV, = 200 MeV, (av) = 
8.2 ■ 10~^^ cm^/s and odm = -04, for an Einasto profile with a = .2 and the med propagation 
parameters. The best fit including PAMELA/ATIC/PPB-BETS without HESS is plotted 
in orange. The parameters are = 900 GeV, = 200 MeV, (av) = 1.1 • 10~^^cm^/s, 
with no parton shower, an Einasto density profile with a = .17 and MIN propagation. 



ters. Additionally, for comparison. Figure 15 shows the same parameters except we change 
the DM profile to Einasto a = 0.17. Furthermore, we show in the figure how HESS al- 
ters this fit by giving the best fit case where only PAMELA and ATIC/PPB-BETS are 
included in the x^- This corresponds to a best fit point of = 900 GeV, = 200 MeV, 
(cTf ) = 1.1 • 10~^^ cm^/s, with no parton shower, an Einasto density profile with a = 0.17 and 
MIN propagation parameters. As we can see, one can find a fit for annihilating background 
that includes a photon background for HESS and satisfies the data. Still, for a good fit, 
the DM profile is required to be less cuspy than the preferred profile of a = 0.17. We can 
also quantify a suppression factor for an Einasto profile with a = 0.17, that could come in 
principle from several sources, so that the effective J is as good as our best fit point. For 
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= 200 MeV one needs a suppression factor of 3-4, for m,^ = 500 MeV we need a factor 
of 4, while for 1.2 GeV we need a larger suppression of 0(10). 

A recent paper [16] investigated the case of DM annihilating through a light and 
then into leptons. Their claim was that models of this type that satisfied PAMELA/ATIC 
with standard DM density profiles were ruled out unless an order of magnitude local boost 
factor was included. Their most stringent bounds come from examining Sgr A*, in radio 
frequencies similar to [15]. Unfortunately, such a constraint is highly sensitive to knowledge 
of the DM density profile at distances much below 100 pc which is beyond the resolution 
limit of current simulations. Nonetheless, if one takes the bound at face value we find that 
their band that agrees with the PAMELA measurement is different than ours, and we can fit 
both PAMELA and their radio bound. In the region of high energy gammas this constraint 
is not more significant than the HESS GC one. The authors also consider bounds from the 
the Sagitarius dwarf galaxy. Since the Sagitarrius dwarf galaxy is being tidally disrupted 
by the Milky Way it is dubious to trust bounds coming from this alone. We have checked 
for one of the most studied dwarf galaxies, Draco (that is not too close to our Galaxy to be 
tidally disrupted), that the bounds are easily satisfied by at least an order of magnitude. As 
discussed, the most significant bound we find arises from the study of the GR region, and 
as shown in Figure [15] they can be satisfied. 

While there exist points in parameter space that can satisfy all experimental constraints, 
this does not mean that generically there is no tension between DM annihilation models 
of this type and the experimental results. As we can see from Section |4| in many regions 
of parameter space the HESS experiment would completely rule out models from the GR 
data. This is due to the fact that the HESS experiment in both the GC and the GR, 
records data that has a best fit to a power law. Once one includes DM annihilating in 
the energy range that HESS studies, it automatically introduces a non-power law shape 
on top of the background. Since the PAMELA experiment currently studies energies less 
than the photon energies recorded by HESS, there is no tension between these experiments. 
However, the inclusion of ATIC/PPB-BETS experiments, automatically signal a mass scale 
which creates the tension. We demonstrate this in Figure 16 where we separate our 
function into Xhess Xatic/pamela P^°^ dimensional contours of for each, 
in the space of and (av) . As we can see from Figure 16 the best fit for taking 
into account PAMELA/ ATIC/PPB-BETS alone prefers a lower mass scale and higher cross- 
sections, while HESS prefers a low cross-sections and potentially higher masses in order to 
preserve its pure power law background. 

This generic tension could mean several things. First, the excesses in PAMELA and 
ATIC/PPB-BETS may simply be caused by astrophysics. Pulsars for instance could ex- 
plain the leptonic excesses without necessarily introducing a large component of high energy 
gamma rays at the center of our galaxy. Second, it could be a red herring, and as we have 
demonstrated, models of this type do have points in parameter space that could account for 
the experimental data. Another way to avoid the tension would be to investigate models of 
decaying DM [7] instead of annihilating DM since the amount of photons at the center of 
the galaxy would then scale like p instead of p^. While these are all possible interpretations. 
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Figure 16: Comparison between the Xhess Xatic/pamela 
of mass and cross-section. 
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we wish to point out another. 

The starting point of current DM investigations is to automatically assume all the recent 
excesses reported by experiments, and work under the assumption that there is/will be no 
conflicting data. For instance in the case of the ATIC/PPB-BETS experiments the excesses 
that are shown conflict with the existing EC [78] data. It has been pointed out in the 
past that EC was based on a much smaller detector area and thus could be more prone 
to systematic errors. This could be the cause of the apparent discrepancy between the 
ATIC/PPB-BETS and EC data. Still, it is curious nonetheless that ATIC/PPB-BETS and 
HESS are the experiments that seem to have a tension if the picture of annihilating DM is 
true. One possibility could be that the anomaly observed by ATIC disappears in the future, 
or that the excess they see is unrelated to DM but rather to astrophysics (e.g. cooling). In 
this case a model would only need to satisfy the PAMELA experiment and the mass scale 
could be much lower in principle. 



In Figure 17 we demonstrate the consequences of assuming only the PAMELA experiment 
and we fit to the electron flux below 100 GeV. The fact that this is a viable option is not 
immediately clear without calculating the results for the other experiments. Indeed, if one 
lowers the mass scale, there are several potential difficulties that need to be confronted. One 
needs to make sure that no feature is introduced in the electron flux that would have been 
seen thus far. At lower energies, even though one can avoid the bounds from the HESS GR 
data, in principle one could now be in danger of over-populating the low energy photons for 
the EGRET experiment. As we show in Figure 17 for the case of ~ 250 GeV, 
MeV and (av) = 9.65 ■ 10~^^cm^/s, with an Einasto density profile for a 
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37 



PAMELA Positrons 



Leplon Flux 



-? 2.0 - 



= 200 MeV 
= 500 MeV 
= 1 ,2 GeV 



i> [GeV] 
HESS Galaclic Ridge 



5x10-^ 



1 X 10-^ 



^ SxlO-' 



; X 10-" 

5x10-' 



50 100 200 

Ey [GeV] 



500 1000 



Up-^ohig Muons from in SK 



Cone Half Aiigie [deg] 



Figure 17: We plot the best fit for = 0.2,0.5,1.2 GeV excluding the data from 
ATIC/PPB-BETS above 100 GeV. The best fit point corresponds to = 250 GeV, 
= 200 MeV and (av) = 9.65 ■ lO^^^cm^/s, with an Einasto density profile for a = 0.17 
and MED propagation. 

propagation, one can satisfy all experimental constraints easily with x^/dof ~ 1. This is 
an incredibly interesting prospect as it opens up possible new avenues for model building 
and additionally provides even more testable predictions. By lowering the mass scale the 
high energy photons are now completely within reach of the FERMI experiment [17] and 
furthermore a nontrivial turnover in the positron fraction measured by PAMELA could be 



observed. In Figure 17 we plot the photon fiux expected from the high energy photons 
observable by FERMI when looking at the GC. This is well within the sensitivity range 
of the experiment, and it is an open question whether or not the high energy photons are 
observable compared to the ICS contribution. We postpone this for future work [79]. 

As we have stated, there are several possibilities for explaining the excesses: annihilating 
DM for low or high m^, decaying DM, Astrophysics/pulsars, some unknown idea or combi- 
nation thereof. Whatever it turns out to be, we are in a golden age of experiments and it's 
useful to review what are the most important experimental results/possible measurements 
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that could be done to sort out potential DM candidates. 

• PAMELA: Future measurements of the positron flux can tell us whether or not the 
fraction continues to increase throughout its mass range and has bearing on whether 
or not DM is heavy or not. 

• FERMI: Has the ability to confirm or reject the ATIC/PPB-BETS excesses by mea- 
suring the electron flux. Additionally if the DM is light FERMI will be in the exact 
range needed to study it's properties. 

• ATIC: The release of ATIC-4 data will offer better statistics and will allow for a 
better comparison to HESS's recent release of the lepton flux at the high end of ATIC's 
reach [80]. 

• HESS: Has already bounded high scale annihilating DM by studying the GR with 
only a relatively short amount of data taking. In principle by collecting more data, it 
could rule out beyond a shadow of the doubt high scale annihilating DM. 

These possible results combined with other experiments should allow us within the next 
few years to confirm or rule out many possibilities for DM. As it stands now, we have 
demonstrated that the particle physics module we have implemented can account for existing 
data. In the near future we will hopefully be able to further pin down or rule out the 
properties of models that have the features we examined. 



Note Added: While this paper was in preparation we learned of a similar work in 
progress by Jeremy Mardon, Yasunori Nomura, Daniel Stolarski, and Jesse Thaler [81]. 
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A The Leaky Box Approximation 

In this Appendix we describe the leaky box approximation which is used to estimate the 
background for the positron flux for energies < 100 GeV. This model is essentially a simplified 



version of that described in section 3.2.1 Its value lies in its minimal set of assumptions, 
allowing one to examine the data with as few extra theoretical inputs as possible. 

The confinement time of charged particles in the Galactic disk is enhanced due to galactic 
magnetic fields and is of order tesc ~ 10^(i?/GeV)~°'^ yr. On the other hand, the cooling 
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time due to ICS and synchrotron radiation is of order tcooi ~ 2 x 10^(i?/GeV)~^ yr. Thus for 
energies below ~ 100 GeV, the particles either reach the Earth or escape the Galaxy before 
loosing energy and therefore cooling can be neglected. For the background computation 
below, we assume this is the case. 

The leaky box approximation is a simplified diffusion model which takes into account 
the confinement of charged particles. The model assumes a free homogeneous diffusion of 
charged particles within the galactic disk. At the galactic boundaries, particles are either 
reflected or escape with finite energy-dependent probability. In its simplest form, the only 
independent parameter is the mean density of matter, Aesc, traversed by the changed particle 
before escaping. One has Aesc = PisuPctesc where pisM is the average interstellar matter (ISM) 
density in the galaxy (not to be confused with the density in the galactic disk) and tcsc is the 
escaping time. Under the assumption that Aesc depends only on the rigidity, 71 = pc/Ze of 
the charged particle, Aesc is extracted by measuring secondary to primary ratios. The B/C 
measurements from the HEAO-3 experiment [48] and other balloon experiments (the most 
recent being ATIC-2 [49]), provide the most stringent constraint. For particles with rigidity 
7^ > 4.5 GV one finds [82,83], 

Aesc = 23.8 P g cm-2. (A.l) 



with S ~ 0.6. As discussed in section 3.2.1 other propagation models allow variations in S 
between 6 ~ .45 — .85. To be conservative we consider these values here. 

In the absence of cooling effects and losses due to collisions, the master transport equation 
at equilibrium for a stable nuclei takes a simple form, 

'^ = Qf'^iE,) + QriE,)- (A.2) 
Here qp™^*^'^^) jg ^i^e primary (secondary) source term for a particle of type i with. 



gsec(^^) ^ ^^PISM ^ f ^^^^^ , daiE, ^ E,) 



(A.3) 



Hi is the number density per unit energy which is related to the flux through $i = (t'j/47r)rij. 



For positrons, where primary sources are absent, the first term on the RHS of (A.2) vanishes 
and one finds, 

n,+ {E) = Q'^{E)tUE) = ^{ripap), (A.4) 

were Q^^+ is the secondary source of positrons generated at the point of interaction, and 
{upOp) = J2j I dEj nj{Ej) -^{Ej — > Ei). For positrons spallation occurs through interac- 
tions of protons and a particles with ISM. Since to a good approximation dp does not depend 
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on energy (other than a sharp kinematical cutoff [84]), a prediction of this theory is that the 
ratio of positron-to-proton flux is proportional to the escape time, 



oc An 



GV 



(A.5) 



The positron injection spectrum, Q^J:, can be computed using the measurements of the 



protons flux which leads to the fit given in Eq. (4.4) and the table below it. The differential 



cross-section for the spallation was derived in [84] and the theoretical uncertainties were 
analyzed in [28]. In light of these uncertainties, we give a rough estimate for the positron 
source (taking for positrons eTZ ~ E), 



Qe 



4 X 10^2^ 



(— V 

VGevy 



-2.84±0.02 



cm ^ s ^ sr ^ GeV \ 



(A.6) 



were we took pism/'^tt. = 0.9 cm ^. Using eq. (A.l), the results, Eq. (3.15) follows. We stress 



that the coefficient for the positron flux may suffer from corrections of order 100%. 
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